On a possible photon origin of the most-energetic AGASA events 
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In this work the ultra high energy cosmic ray events recorded by the AGASA experiment are analysed. With 
detailed simulations of the extensive air showers initiated by photons, the probabilities are determined of the 
photonic origin of the 6 AGASA events for which the muon densities were measured and the reconstructed 
energies exceeded 10 20 eV. On this basis a new, preliminary upper limit on the photon fraction in cosmic rays 
above 10 20 eV is derived and compared to the predictions of exemplary top-down cosmic-ray origin models. 



1. Photons as UHECR? 

Studies on ultra-high energy photons as pri- 
mary cosmic rays have the potential to discern be- 
tween different models of origin of ultra-high en- 
ergy cosmic rays (UHECR) . Many top-down sce- 
narios predict a dominance of photons in the CR 
flux above 10 20 eV (for a review see while ac- 
cording to classical acceleration scenarios in this 
energy range the CR flux should consist mainly 
of hadrons. Thus, any experimental conclusions 
about a fraction of photons in UHECR, e.g. an 
upper limit on this fraction, are expected to pro- 
vide constraints for the models of UHECR origin. 

Presently, there are few experimental upper 
limits on photon fraction estimated on the ba- 
sis of cosmic rays of energies above 10 19 eV (£3], 
They are plotted in Fig. ^ an d compared to 
the predictions of the three exemplary top-down 
models: super heavy dark matter (SHDM) 0], 
Z-bursts (ZB) and topological defects (TD) 0. 
It can be seen that the available limits are not 
very constraining since below 5xl0 19 eV a typi- 
cal predicted photon fraction is still low. Above 
10 20 eV, where the photon fraction predictions 
are larger, there is no experimental upper limit, 
although there are a few well documented cos- 
mic ray events of such extremely large energies. 
While large amount of UHECR data from next 
generation giant cosmic-ray detectors is still on 
their way, in this work the available data of ex- 



tensive air showers detected by the AGASA ex- 
periment are analysed to derive an upper limit 
on the photon fraction in UHECR. In order to 
determine the probability of an event being ini- 
tiated by a photon the muon densities measured 
on ground, which are expected to be good indica- 
tors of primary mass |6I7| . were compared to the 
simulations. 

2. AGASA cosmic ray data above 10 20 eV 

Among the 11 cosmic ray events above 10 20 eV 
recorded by the AGASA experiment |8I2| in 6 
cases the muon densities at a distance of 1000 m 
from the shower core (p M (1000)) were measured 
and published in Ref. (see Table The accu- 
racy of p^(1000) was evaluated to be ~ 40% 

2.1. Simulations 

The published energies and arrival directions 
were used for the simulation of photon-induced 
air showers. The effect of pre-cascading of pho- 
tons in the geomagnetic field (preshower effect) 
1101 and the subsequent air shower simulations 
were performed with use of the CORSIKA pack- 
age Q3] with coupled PRESHOWER code [E]. 
For each of the 6 events with available experimen- 
tal ^(1000) values, 100 photon-induced showers 
were simulated. This allowed obtaining 6 distri- 
butions of simulated muon densities p M . S im(1000) 
at the observational level. 

The exact values of the reconstructed primary 
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energies E were published in Ref. ^3] an d up- 
dated on the AGAS A web page |21 ■ These values 
were obtained with accuracy of 25%, assuming 
that all the showers were initiated by hadrons. 
As noted in Ref. 2 , for photon-induced show- 
ers, Eq underestimates the real primary energy 
by about 20%. This 20% correction was applied 
in all the photon-induced shower simulations pre- 
sented here. 

2.2. Probabilities of photonic origin 

The basic approach to determine the prob- 
ability of a primary photon origin of an indi- 
vidual event is quite similar to the one out- 
lined in Ref. ^Hl- The simulated distributions of 
P/j, sim (1000) were compared to the experimental 
values in the following way. It is assumed that 
the uncertainties connected to the measurements 
of the muon density ,0^(1000) for each AGASA 
event follow the Gaussian distribution with the 
uncertainty <t P(i (iooo) = 0.4/3^(1000). The simula- 
tions for the energy and arrival direction of the 
j-th recorded event gave n — 100 expected values 
/9^ sim (1000), where index i denotes the i-th sim- 
ulated shower. The probability P lJ that /^(lOOO) 
differs by more than S lJ from the expected muon 
density /Q^ S i m (1000) is given by 

P ij = P(|p£»m(1000) - ^(1000)| > «F) = 

[Pi (1000)+<5 ij ' 

= 1- f(x,pl(1000))dx, (1) 

with f(x, p^(1000)) being the Gaussian distribu- 
tion. Since cr Pfi rs 0.4p M , / has a small (less than 
1% of J f(x)dx) contribution from negative val- 
ues of muon densities x. This contribution was 
accounted for by a simple renormalization of P iJ . 
The probability for the j-th event being initiated 
by an UHE photon, P$ was obtained by taking 
the average P 4 -? values. The experimental data 
together with the results of the simulations and 
analysis are collected in Tabled F° r each event, 
the probabilities P^ are small, but not negligible. 

3. Upper limit on photon fraction in 
UHECR above 10 20 EeV 

After the probabilities p? have been deter- 
mined, the upper limit for the photon fraction 



Table 1 

Average muon densities (p^ sim (1000)) and their 
RMS fluctuations simulated at ground level for 
the AGASA events with E > 10 20 eV are com- 
pared to the experimental values p^(1000). The 
Eq values are increased by 20% with respect to 
the published data. The probabilities P^ that 
the j-th event was initiated by an UHE photon 
are given. The azimuths of the arrival directions 
increase clockwise from North. For the analysis 
method see the text. 

AGASA data 

E [EeV] 295 240 173 161 126 125 

zenith [°] 37 23 14 35 33 37 

azimuth [°] 260 236 211 55 108 279 

^(1000) [m- 2 ] 8.9 10.7 8.7 5.9 12.6 9.3 
simulations 

preshowers [%] 100 100 96 100 93 100 

pj> s (1000) [m~ 2 ] 4.3 3.1 2.1 2.3 1.7 1.8 

RMS(p£ s )[m- 2 ] 0.7 0.6 0.8 0.3 0.4 0.3 

P? [%] 19.5 7.3 6.0 12.3 2.4 3.8 

Table 2 

Probability that A 7 out of the 6 AGASA events 

were initiated by a photon. 
A 7 1 2 3 4 5 6 
P[%] 57.8 34.0 7.4 0.8 4E-2 1E-3 1E-5 

in the cosmic rays with Eo > 100 EeV can be 
estimated in the following way. Let us calculate 
the probability that 7V 7 out of 6 AGASA events 
were induced by photons: 

n 6— n 

P(A 7 =n)=^n^ fc Il( 1 - P 7) ( 2 ) 

fe=l 1=1 

where j k> ii € A = {1, 6}, V k, I G A : j k ^ ji 
and the summation extends over all n-elemental 
combinations of indices j k in A. The values of 
P(A 7 = n) are collected in Table |21 These val- 
ues determine the upper limit on gamma frac- 
tion in cosmic rays. The probability that less 
than 3 out of the analysed 6 AGASA events with 
Eq > 100 EeV were initiated by photons is larger 
than 95%, i.e. P{N 1 /N < 3/6) > 95%. This al- 
lows the conclusion that the upper limit on the 
photon fraction in cosmic rays above 10 20 eV 
equals 3/6 = 50% at 95% confidence level. 
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Figure 1. The ratio of photon- induced showers to 
hadron-induced ones as a function of the primary 
particle energy 2 . The upper limits obtained by 
the AGASA and Haverah Park experiments at 
95% CL are indicated by arrows. For explanation 
of the theoretical models see the text. 

4. Discussion 

The obtained 50% upper limit on gamma frac- 
tion corresponds to N 1 /N nuc i(> 10 20 eV) < 1. In 
Fig. Q the above result is compared to the theo- 
retical expectations of exemplary top-down mod- 
els of cosmic ray origin and to the previous ex- 
perimental upper limits on the photon fraction in 
UHECR of energies above 10 20 eV. The models 
that predict gamma fractions higher than the ob- 
tained upper limit are disfavored by the AGASA 
data. 

It should be noted that a considerable system- 
atic uncertainty is related to extrapolating the 
photonuclear cross-section er 7 _ Q i r to highest en- 
ergies |16|. In relation with the current analysis, 
assuming that a 7 - a i r increases more rapidly at 
highest energies would result in larger muon den- 
sities on ground and hence in larger probabilities 
of photonic origin of the analysed events. This 
will be discussed elsewhere. 

The presented analysis method of determining 
the upper limit of photon fraction can be used as 
well for a larger data sample which is expected 
soon from the new generation cosmic ray detec- 
tors. The conclusions based on these new data 
should offer a good way to distinguish between 
different models of cosmic ray origin. 
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